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Ci POTENTIAL HARDWARE AND SOFTWARE IMPROVEMENTS OF 

ZD INERTIAL POSITIONING AND GRAVITY VECTOR DETERMINATION* 


O 

rH 

< 



H, BAUSSUS VON LUETZOW 


ZUSAMMENFASSUNG 

In dlesem Vortrag werden der potentlelle Elnschluss von Krelseln und 
Beschleunlgungsmessern hoeherer Genauigkelt in den horlzontalen Kanaelen 
des vertikal orientlerten Inertialortsbestimmungssystems von Litton 
und resultierende Verbesserungen hlnslchtllch Posltlons-und Lotabwelch- 
ungsennittlung eroertert. Ferner werden Lotabweichungs-Fehlerschaetz- 
ungen mittels einfacher und hoeherer Datenreduzierungsverfahren 
geliefert und interessferende Beschraenkungen Identifiziert. 


SUMMARY 

The paper discusses the potential incorporation of higher performance 
gyroscopes and accelerometers in the horizontal channels of Litton's 
local-level Inertial positioning system and the resulting improvements 
in positioning and deflection of the vertical determination. It 
further provides deflection error estimates based on the use of simple 
and advanced data reduction methods and Identifies limitations of 
interest. 

RESUME 


Cette etude est une discussion sur 1'incorporation potentlelle des 
gyroscopes de haut rendement et des accelerometres dans les voles 
horizontalles du systeme Inertlelle de determination du point a niveau 
locale de la compagnle Litton. Cette communication fournit une 
evaluation de l'erreur de deflection basees sur l'emploi de m§thodes 
fondamentalles et avancSes de reduction de donnees. Les limites de 
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TRODUCTION 

e present Rapid Geodetic Survey System (RGSS), developed by Litton 
stems, Inc. for the U.S. Army Engineer Topographic Laboratories (ETL), 
s the capability of determining horizontal and vertical positions, 
flections of the vertical, and gravity anomalies with average rms 
rors of 1 m, 0.3m, Earcsec, and 2 mgal, respectively, for 50 km 
ns under utilization of post-mission adjustments. It operates as 
quasi local-level system. It thus does.not require altitude 
mplng, permits Kalman stochastic enrol* control without great 
mplexlty under consideration of observed velocity errors at 
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vehicle stops, .and employs effective post-mission adjustments 
with the aid of terminal position and gravity vector Information 
only. Present critical hardware consists of an A-1000 vertical 
accelerometer, two A-200 horizontal accelerometers, and two 
6-300 gyroscopes. Although data utilization from repeated runs 
has resulted In average deflection rms errors of about 1 arcsec, 
the requirement of a maximum rms error of 0.5 arcsec necessitates 
the Incorporation of higher performance gyroscopes and acceler¬ 
ometers. This would simultaneously achieve conventional surveying 
accuracy (10~ 5 ). 

2. SIGNIFICANCE OF HIGHER PERFORMANCE GYROSCOPES AND 
ACCELEROMETERS FOR POSITION AND VERTICAL DEFLECTION 
ERROR REDUCTION 

Following the approximate elimination of the effects of constant 
gyro biases on positions under utilization of accurate Initial 
and terminal coordinates, the remaining dominant error source 
is gyro correlated random noise. In this respect, the pertinent 
parameters for the 6-300 instruments are a standard deviation 
of 0.002° hr -1 and a correlation time of 3 hours. Another error 
source to be considered is the accelerometer scale factor, 
assumed to be constant for a test run. The standard scale 
factor error for the A-200 accelerometer Is 0.01%. Correlated 
accelerometer noise Is characterized by a standard deviation of 
10 mgal and a correlation time of 40 minutes. Both correlated gyro 
random drift and correlated accelerometer measurement errors 
affect the accuracy of deflections of the vertical estimated 
from Inertial data and Initial and terminal deflection components, 
possibly augmented by corresponding geodetic azimuths. In 
order to achieve deflection accuracies smaller than or equal 
to 0.5 arcsec rms. It is, therefore, necessary to Install higher 
performance gyros and accelerometers. A significant reduction 
of the two autocorrelation parameters is also a necessary 
prerequisite for the optimal statistical estimation of deflections 
of the vertical. The Litton 6-1200 gyro Is expected to be 
compatible with the stringent deflection accuracy requirement. 

Its random drift error is about O.OOl^r- 1 and its autocorrelation 
parameter Is short. For details, reference Is made to 
Litton [1974] and Huddle [1977a].The A-1000 accelerometer 
addressed by Litton [1973,1975] has a standard scale factor error 
of approximately 0.005%. Correlated accelerometer noise Is 
2 mgai rms, equivalent to aboutO.35 arcsec. As mentioned In 
the Introduction, this Instrument Is already Incorporated In the 
vertical channel of the present RGSS. In conjunction herewith, 
an Improved velocity quantizer has been Installed. These 
Improvements have resulted in altitude errors of 0.3m rms 
and gravity errors of 2 mgal rms after 2-hour runs and under 
utilization of post-mission adjustments. Of further 
significance,the precision of the 6-300 gyros and A-200 
accelerometers, primary horizontal channel components of the 
present R6SS, Is greater than Indicated by the accuracy 
parameters. This has been observed by Schwarz [1979*1 and 
by ETL after evaluation of repeated test runs during 1980. 
Multiple runs can also be expected to reduce the single run 
standard deviations associated with an advanced RGSS (AR6SS), 
advocated by the author in 1978. * 
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3. BASIC METHOD FOR THE DETERMINATION OF DEFLECTIONS 
OF THE VERTICAL 

A presentation of both Kalman error control and post-mission 
adjustment and of the determination of deflection components has 
been made by Baussus von Luetzow [1979], With respect to position 
determination, Schwarz [1979b] noted that the errors Increased 
considerably during tests In strongly mountainous terrain. This 
Is due to the fact that under these conditions gravity anomalies 
and deflections of the vertical cannot be satisfactorily modeled 
as quasl-statlonary random variables In the real time Kalman 
error controller. Fortunately, these Inherent difficulties do 
not significantly affect the determination of deflections of 
the vertical,the detailed discussion of which requires the 
utilization of the local-level error differential equations, 
applicable to horizontal motion and restricted to land vehicles: 
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+ w z»k + Vz + y 
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For simplicity, the symbol 6 In front of the dependent variables 
has been ommltted. The applicable coordinate system Is evident 
from figure 1. 

x (north) 


A 


> y (east) 


Figure 1 

Applicable Coordinate System 
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Symbols used In the foregoing equations are with z downward: 

x north angular position error 

y east angular position error (for longitudes west of 
Greenwich) 

x north velocity error 

y east velocity error 

a azimuth axis angular drift rate error 

B north axis angular drift rate error 

y east axis angular drift rate error 

g normal gravity 

f geographic latitude 

R mean earth radius 

® z azimuth platform attitude error 

» N platform tilt error about north axis 

*£ platform tilt error about east axis 

gn product of g and deflection component n 

g£ product of g and deflection component & 

dVv 

V—a north acceleration of survey vehicle 

dV F 

S «—£ east acceleration of survey vehicle 
E dt 

a E correlated east accelerometer error 

a N correlated north accelerometer error 

n N - ncos4> north earth rate 

p n .V e /R north angular rate 

“N* n N + p m ■ Hcos4 + V e /R north spatial rate 

«£■ p E * -V„/R east spatial rate 

+ p 2 » flsln* 4 v e /R • tn$ vertical spatial rate 

In Inertial land navigation terms involving » E in equations (5) 

and (6) are neglected In Lltton's Kalman error controller. The 
Initial conditions at time to ■ 0 are 
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In the general RGSS node of operation, position and gravity vector 
Information Is generated. The corresponding Kalman mechanization 
Is then simplified by setting ■ u» * 0 In equations (5)-(7) 
which results In the uncoupling of tne horizontal channel error 
differential equations. This Is not critical for travel periods 
of 2 hours. 

As an example, equations (4) and (7) may then be formulated at 
the end of the first travel period, whence t ■ ti. 


* 



- *i ♦ a Ni 

- X I + f vdt ♦ e 0 
r J 




( 8 ) 

(9) 


In eq. (8), the existence of an Initial accelerometer measure¬ 
ment error has been taken Into consideration. 

The magnitude Xi Is well estimated by the observed velocity error 
xi. After Implementation of the tilt correction (est.), equations 
(9) and (8) read 

•i}' 1 “ -«i + J \ dt + C 0 (10) 


* 1 ° « g(-«i - g^i - C 0 \ + a Nl - a N() (11) 


The estimation error Is generally negligible. 

Because of the quasl-ellmlnatlon of Integrated terms Involving 
A, y% x In equations (5)-(7) by appropriate tilt corrections, 
constant gyro biases may be obtained with a high degree of „ 
approximation by a solution of these equations with x » y « x * 0 
and o ■ 5, 6 ■ I* y • Y» It Is then, according to Huddle [1977bJ, 

*E " *(* + ® Y*) ( 12 ) 

£* t + B ^ (13) 

'■■'('T* * ‘ ) 


( 14 ) 
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Under consideration of terminal deflection and azimuth closure 
errors. It Is possible to compute a, 8, 7 from 


* * estimated " * observed 


n estimated * n observed 


^ estimated ^ observed (^) 


where A designates geodetic azimuth. 

Examination of equations (12)-(14) reveals that In moderate 
latitudes there exists relatively small coupling for time Intervals 
not exceeding 2 hours. For this reason, linear approximations 
♦e«* yt and 8t have been used with success. Highly accurate 
deflection determinations, accomplished by means of an ARGSS, would, 
however, require consideration of a terminal azimuth error. 

For t * tj, equation (11) may be formulated with * 0 as 

x(ti) « x ± - g? El - 9 Ui - i 0 ) + g / 1 (v-7)dt + a Ni - (18) 


i(t n ) - x n « 9®^ - 9 U n - 6 e ) ♦ 9 / (y-T)dt + #Nn - a No (19) 

o 


Equations (12) and (19) show that the computation of * Ei is to 
the first order associated with an error 

5? Ei * - r C 9 / tn(Y ‘ 7)dt + a Nn * a N 0 3 (20) 

XI o 


The basic solution for the prime deflection of the vertical from 
equations (18) and (19) Is 



( 21 ) 



Equation (21) may be supplemented by 




«o + £ 


(6C- - - (*. * 7 1 *e ) 

e ° s i tn Sn 


( 22 ) 


to account for astrogeodetlc deflection errors and accelerometer 
bias errors. For a straight traverse, the last term In eq. (22) 
tends to cancel out. 

The rms deflection error orfc^) can be computed by covariance 
analysis Involving the terms without parentheses in eq. (21). 

Under Inclusion of the first two terms of eq. (22) It Is 

var * var 4 var y* + 0 - |^) 2 var £ 0 4 (51) 2 var (23) 


where var ai Is the accelerometer-Induced variance and var 
designates the gyro-induced variance. Under consideration of present 
RGSS parameters. Identified in section 2, the rms deflection error 
c£. is approximately presented In Figure 2. 



Distance Traveled (km) 


Figure 2 

Approximate RMS Deflection Error as Function of Travel Time 


The variation of normal gravity g with altitude Is not critical 
In terrestrial applications. The scale factor-generated error 
does not vanish with respect to L-shaped traverses. It Is, 
however, relatively small for travel times not exceeding 2 hours. 
Heading sensitivity Induced by significant azimuth changes 
appears to be somewhat critical and may require empirical 
corrections. Schwarz [1979b] discussed the Impact of heading 
sensitivity on positioning and concluded that error reductions may 
be achieved by a modified smoothing procedure. The problem Is more 



610.1/8 


Intricate In the case of deflection determinations. The use of 
static accelerometer measurements In conjunction with Initial 
and terminal deflections together with a simplified Kalman 
filter originated with Huddle [1977a3. The above analytical 
presentation Is one of several quasl-llnearizatlon techniques 
according to the terminology of Schwarz [19803. It may be 
refined by a finite difference solution of the system 
differential equations (l)-(7) with + • const, with or without 
periodic filter corrections. This would result In weight factor 
expressions as a function of time involving the generating 
variables e • C[x(t), y(t)3, n • n[x(t), y(t)3, » E ' « E • ~ *n • 

a, 6, and y. Approximate constant vehicle velocities would f 
simplify the computations. A matrix solution involving all 
unknown c s and n’s would then replace the independent single 
component solutions. 

4. ADVANCED METHOD FOR DEFLECTION DETERMINATION 

Equation (21) under Inclusion of the first two terms of eq. (22) 
or a comparable refined solution Identified In section 3 may be 
formulated as 

Ej * i t * \ (24) 

mm 

where E i * + T £l - ^ Is a message variable, 

Is a signal variable, and the remaining terms denote a noise variable - n lt 
The collocation method in physical geodesy permits In semi-flat 
terrain the estimation 

E e * E «i(Ei + n t ) « Ajfij 4 N t ) (25) 

where Is the matrix of regression coefficients a t to be computed. 

It Is then in matrix form, with bars Indicating covariances, 

+ *k> - 0, 1, ...n (26) 


The solution for the regression coefficient matrix follows as 

♦ W’ 1 < 27 > 


Due to the sizable gyro and accelerometer correlation times associated 
with the present RGS5, the advanced method characterized by equations 
(25) and (27) does not provide significantly better deflection 
estimates. Considerably Improved estimates would, however, result 
under utilization of Intermediate deflection constraints. In the 
context of an ARGSS, l.e., under consideration of very short correla¬ 
tion times, the advanced method may be expected to yield better 
estimates. In addition. It constitutes the statistical framework for 
an optimal area solution under utilization of data relating to 
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several traverses. The computation of topographic corrections and 
of spatial covariance functions for applications In strongly 
mountainous terrain has been presented by Baussus von Luetzow [1981], 
Alternatively, It would be necessary to conduct repeated runs. In 
this way, rms deflection errors not to exceed 0.35 arcsec are within 
the reach of an ARGSS. 

5. CONCLUSION 

The Incorporation of Identified or other higher performance gyros and 
accelerometers and of Improved velocity quantizers in Litton*s RGSS 
would result in an ARGSS. Non-linear gyro bias corrections, 
elimination of errors caused by gyro heading sensitivity, and 
advanced data reduction techniques or repeated runs are then expected 
to facilitate deflection determinations with standard errors not 
to exceed 0.35 arcsec for a travel time of 2 hours or about 60 km 
distance. As a by-product, conventional surveying accuracy would 
be achieved. Deflection changes £-£ 0 and n-n 0 with standard errors 
of the order 0.1 arcsec for the same travel time and vehicle speed 
would require supplementation by a high-performance gravity 
gradiometer. 
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